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APPLICATION  OF  S1BYL2  TO  THE  AGARD  WG18  COMPRESSOR  TEST  CASES 

by 

W.  J.  Calvert 

SUMMARY 

S1BYL2  is  an  inviscid-viscous  blade-to-blade  method  for  calculating  the 
detailed  aerodynamics  and  overall  performance  of  compressor  blades.  It  may  be 
applied  either  on  its  own  to  predict  the  flow  for  individual  blade  sections,  such 
as  the  mid  span  of  a  linear  cascade,  or  in  conjunction  with  a  throughflow  calcu¬ 
lation  to  predict  the  performance  of  a  complete  axial  compressor. 

A  previous  AGARD  paper  by  the  author  described  applications  of  S1BYL2  to 
most  of  the  compressor  cascade  test  cases  which  have  subsequently  been  selected 
by  AGARD  Working  Group  1 8 :•  generally  good  agreement  was  obtained.  This  current 
paper  presents  new  predictions  for  the  V2  and  ARL  SLI9  cascades  and  for  the  high 
speed  compressor  cases.  It  is  hoped  that  this  will  be  one  of  many  sets  of  calcu¬ 
lations  for  these  cases,  so  that  an  improved  understanding  of  each  case  may  be 
obtained,  together  with  an  appreciation  of  the  strengths  and  weaknesses  of 
different  computational  approaches. 
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Application  of  S1BYL2  to  the  AGARD  W31 8  Compressor  Test  Cases 

W  J  Calvert 
Propulsion  Department 
Royal  Aerospace  Establirhrr.^iit 
Pyestock,  Famborough, 

Hampshire,  UK 


SIMiARY 

1BYL2  is  an  inviscid-viscous  blade-to-blade 
u itnod  for  calculat.  vj  the  detailed  aero¬ 
dynamics  and  overall  performance  of  compres¬ 
sor  blades.  X*  may  r>e  applied  either  on  it3 
own  to  predict  trie  flow  for  individual  blade 
sections,  such  as  the  mid  span  of  a  linear 
.as  ade,  or  in  oonji.  !Ction  with  a  throughflow 
calculation  to  predict  th  performance  of  a 
complete  axial  ccnpressor. 

A  previous  AGARD  paper  by  the  author  descri¬ 
bed  applications  of  SIBYLr  to  most  of  the 
compressor  cascade  test  cases  whio.  havr. 
subsequently  been  selected  ry  AGARD  Working 
Group  18  :•  generally  good  agreement  was 
obto.ned.  This  current  paper  presents  new 
pro  li  -t  i  ons  for  the  V.  and  AA1.  SL19  cascades 
and  tor  tne  high  speed  ccnpressor  cases.  It 
is  hovel  that  this  will  be  one  of  many  sets 
of  calculations  for  these  cases,  so  that  an 
improved  understanding  of  each  case  may  be 
obtained,  together  with  an  appreciation  of 
the  strengths  and  weaknesses  of  different 
computational  approaches. 

LIST  OF  SYMBOLS 

c  blade  chord 

H  boundary  layer  shape  factor  8*/0 
M  Mach  number 
P  static  pressure 
a  flow  angle,  degrees 
8*  boundary  layer  displacement  thickness 

0  boundary  layer  momentum  thickness 

0)  loss  coefficient  based  on  upstream 
conditions 

fl  stream  tube  contraction  (upstream/ 
downstream) 

Subscripts 

1  conditions  at  upstream  boundary  of 
cascade  or  inlet  to  blade  row 

1  conditions  at  downstream  boundary  of 
cascade  or  exit  frctn  blade  row 

1.  INTRODUCTION 

Validation/calibration  of  computer  codes 
against  suitable  te  ;t  cases  is  a  highly 
important  area  of  Computational  Fluid 
Dynamics  which  is  rightly  receiving  increased 
attention,  as  shown  for  example  by  Ref  1. 

The  current  state-of-the-art  for  turbo¬ 
machinery  test  cases  is  indicated  by  the 
selection  made  by  AGARD  Working  Group  18 
(W318)!.  These  cover  a  useful  range  of 
practical  cascades  and  rotating  machines, 
and  provide  a  reasonable  amount  of  data  on 
the  flow  fields  as  well  as  measurements  of 
overall  performance.  However,  as  noted  by 
W318  in  Chapter  III,  neither  the  experi¬ 
mental  data  nor  the  flow  codas  are  perfect. 
Also  it  is  not  generally  possible  to  obtain 


experimental  data  of  many  important  areas  of 
high  speed  flows,  due  primarily  to  the 
relatively  small  scale  of  turbcrrachinery 
blading.  Thus  sare  aspects  of  the  code3 
cannot  be  adequately  validated.  Given  this 
situation,  there  is  ruch  to  be  gained  by 
applying  codes  to  a  wide  range  of  test  cases 
aid  by  comparing  the  results  iron  different 
codes  for  the  same  test  oases.  An  important 
contribution  of  W318  is  to  have  identified  a 
suitable  set  ol  test  cases  to  initiate  this. 

The  main  code  employed  in  the  present  work  is 
ti.e  RAE  S13TL2  method  3.  This  is  an  invrscid- 
viscous  interaction  technique  for  predicting 
the  blade-to-blade  performance  of  axial 
compressors.  It  can  be  applied  either  on  it3 
own  to  individual  blade  sections,  such  as  the 
mid  span  of  *.  linear  cascade,  or  linked  to  a 
streamline  curvature  throughflow  method  in  an 
S1-S2  system  to  give  a  complete  qua3i-three- 
dimensional  prediction  of  the  flow  and 
overall  performance  of  one  or  more  blade 
rows. 

A  brief  description  of  the  S1BYL2  and  S1-S2 
methods  is  given  in  section  2  of  this  paper. 
Section  3  then  reviews  the  results  from 
previous  applications  »  of  S1BYL2  to  the 
compressor  cascade  test  cases  chosen  by  W318 
and  presents  further  predictions  for  the  V2 
high  subsonic  cascade  and  new  predictions  for 
the  ARL  SL19  transonic  cascade.  Results  from 
the  RAE  S1-S2  system  for  the  tnree  high  speed 
compressor  cases  are  then  given  in  section  4. 

2 .  GENERAL  DESCRIPTION  OF  METHODS 

The  RAE  S1BYL2  method  is  an  inviscid-viscous 
interaction  technique  to  predict  the  blade- 
to-blade  performance  of  axial  compressor.,. 

The  inviscid  part  consists  of  a  time  marching 
Euler  calculation,  based  on  Ref  5,  30  it  can 
handle  transonic  flows:  3hock  waves  can  be 
captured  and  the  losses  due  to  them  can  be 
predicted.  The  calculation  takes  place  on  a 
specified  axisyimetric  stream  surface,  and  it 
includes  the  effects  of  rotation  and  of 
varying  radius  and  stream  tube  thickness  in 
the  axial  direction.  The  viscous  calculation 
is  an  integral  technique  consisting  of  three 
parts  to  estimate  laminar  boundary  layer 
development,  transition  point  and  turbulent 
boundary  layer  development  respectively. 

The  inviscid-viscous  matching  procedure 
employs  mixed  modes  to  allow  a  valid  solution 
to  be  obtained  even  when  there  are  regions  of 
separated  flow  due  to  shock  wave/boundary 
layer  interactions  or  to  excessive  diffusion. 
The  viscous  calculation  continue  in  the  wake 
downstream  of  the  bladed  region,  and  a 
compressible  flow  mixing  calculation  is 
carried  out  on  the  downstream  boundary  to 

Hoforminp  IHSSn  6Xi£  CQIKlitlionS .  ThliS 

method  provides  both  details  of  the  internal 
flows,  including  boundary  layer  development, 
and  predictions  of  the  blade  section  overall 
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performance,  such  as  exit  flow  angie  and 
losses . 

In  order  to  employ  the  S1BYL2  nrethod  to 
predict  the  performance  of  complete  blade 
rows,  It  must  be  linked  with  a  throughflow 
calculation  which  determines  the  spanwise 
variations  of  the  flow  and  ensures  tliat 
radial  equilibrium  is  satisfied.  A  semi¬ 
automatic  S1-S2  procedure  for  this  is  descri¬ 
bed  in  Ref  6.  This  procedure  has  now  been 
automated,  with  both  the  S1BYL2  calculation 
for  the  blade-to-blade  sections  ana  tne 
streamline  curvature  calculation  for  Use  hub- 
to-tip  flow  incorporated  into  one  program 
together  with  the  blade  geometry  ana  data 
link  routines.  The  streamline  curvature 
routine  1  is  an  improved  modular  version  of 
that  previously  employed.  It  allows  curved 
calculating  planes  to  be  soecified  to  give  an 
accurate  match  to  the  leading  and  trailing 
edges  of  each  blade  row.  Spanwise  misting 
terms  have  been  added  to  seme  versions  of  the 
routine,  but  these  were  not  Included  in  the 
present  S1-S2  calculations. 

The  basic  principles  of  the  S1-S2  interaction 
remain  the  same.  Essentially  there  are  three 
parts  to  each  loop:- 

i  the  performance  of  the  blade-to- 
blade  sections  is  predicted  by  S1BYL2  - 
this  involves  calculations  for  up  to  9 
sections  on  each  blade  row. 

ii  the  hub-to-tip  flow  throughout  the 
whole  compressor  is  calculated  by  the  S2 
streamline  curvature  method,  using  the 
pitchwise  averaged  data  from  the  Si 
solutions . 

iii  the  relative  flow  angles  and 
stream  tube  geometry  from  the  SI  and  S2 
parts  of  the  calculation  are  compared  to 
see  whether  the  solution  has  converged  - 
if  not  new  blade  geometry  for  each 
section  is  found  by  calculating  the 
intersections  between  the  stream 
surfaces  frem  the  current  S2  solution 
and  the  blade  aerofoils. 

Target  convergence  criteria  for  the  complete 
S1-S2  solution  are  typically  a)  SI  and  S2 
relative  flow  angles  agreeing  to  within  0.1°, 
and  b)  SI  and  S2  stream  tube  thickness 
variations  agreeing  to  0.1%.  The  number  of 
S1-S2  loops  required  to  achieve  this  depends 
on  the  difficulty  of  the  blade  section 
operating  conditions  ati  the  number  of  blade 
rows  involved.  Ado.,'-  eug,..,  loops  are 
generally  sufficient  ft.  one  or  two  blade 
rows,  but  more  vail  be  vended,  for  exarrple, 
if  blade  choking  i-  -  qntficantly  affecting 
the  flow  distrir*  cion.  The  solutions  from 
the  previous  loop  provide  the  starting 
conditions  for  calculations  on  the  current 
lap,  so  that  the  procedure  is  quite 
efficient.  For  example,  a  solution  for  the 
Hr 3  cy ik  RL  stage  rsee  section  4.2)  requires 
about  6  hours  CPU  on  a  Stardent  1500  mini- 
supercaiputer . 

It  should  be  noted  that,  in  order  to  obtain 
realistic  losses  towards  the  blade  ends, 
empirically  derived  extra  losses  must  be 


specified  as  an  input  to  the  S1-S2  system. 

These  additional  losses  are  modelled  by  a 
constant  drag  force  in  the  SI  calculation  and 
id  e  effects  are  fed  to  S2  via  the  normal 
S1-S2  linkirg.  Although  the  need  for  such 
input  is  one  of  Uie  limitations  of  the  S1-S2 
approach,  the  ease  with  which  3uch  correc¬ 
tions  can  be  incorporated  in  order  to  tune 
the  modal  for  a  given  compressor  can  be  a 
distinct  advantage  in  practice.  It  i3 
emphasised,  however,  that  apart  from  the 
stated  end  losses  and  annulus  wall  blockage 
factors  no  other  corrections  were  employed 
for  the  cases  presented  here. 

3 .  COMPRESSOR  CASCADE  TEST  CASES 

AGARD  Working  Group  18  selected  four 
compressor  cascades  suitable  for  inviscid- 
viscous  blade-to-blade  methods 

i  E/CA-2  V2  High  Subsonic  Cascade 

ii  E/CA-3  CMERA  115  High  Subsonic 

Cascade 

iii  E/CA-4  EFVLR  L030-4  Low  Supersonic 

Cascade 

iv  E/CA-5  ARL  SL19  Transonic  Cascade 

S1BYL2  ha3  previously  been  applied  to  the 
first  three  of  these  and  the  results  are 
given  in  Ref  4.  There  was  generally  good 
agreement  between  the  test  results  and 
predictions,  in  terms  of  both  internal  flow 
details  such  as  blade  surface  pressure 
distributions  and  overall  performance  para¬ 
meters  such  as  static  pressure  ratio  and  exit 
flow  angle.  The  only  significant  problem 
area  was  that  increases  in  losses  at  high 
incidence  were  not  predicted  directly,  but 
relied  on  the  user  noticing  the  high  effec¬ 
tive  incidence  indicated  by  the  solution  and 
specifying  revised  starting  conditions  for 
the  surface  boundary  layer  calculation.  For 
cases  near  optimum  incidence,  deviation 
angles  were  predicted  to  about  ±  1°,  static 
pressure  ratios  to  ±  4%  and  overall  loss 
coefficients  to  ±  2%  of  inlet  dynamic  head. 

These  levels  of  accuracy  have  been  confirmed 
in  applications  of  tire  code  to  other  cas¬ 
cades,  such  as  the  supercritical  designs 
described  in  Refs  8  and  9. 

The  test  conditions  chosen  by  W318  for  the 
CMERA  115  and  DFVLR  L030-4  cascades  are  the 
same  as,  or  similar  to,  those  previously 
studied.  Therefore  only  calculations  for  the 
V2  and  ARL  SL19  cascades  will  be  presented 
here. 

3.1  V2  Subsonic  Compressor  Cascade 

The  V2  cascade  is  a  high  canter  (56.8°),  high 
solidity  (2.22)  design  with  7%  thick  DCA 
blade  profiles.  Twelve  test  conditions  were 
selected  by  WG18,  covering  inlet  Mach  numbers 
from  0.3  to  0.85,  inlet  flow  angles  from 
47.5°  to  54.5°  and  stream  tube  contraction 
ratios  from  1.08  to  1.39. 

As  noted  in  Refs  2  and  4,  the  axial  variation 
assumed  for  stream  tube  thickness  between  the 
measured  upstream  and  downstream  values  is 
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critical,  but  also  open  to  uncertainty.  Two 
plausibl*.  assumptions  could  be  (a)  a  linear 
variation  with  streamwise  distance  between 
the  two  measuring  planes,  and  (b)  constant 
values  upstream  and  downstream  and  a  linear 
variation  between  the  leading  and  trailing 
edge  planes  of  the  caacacte.  For  the  V2 
cascade,  wnich  has  a  high  aspect  ratio  of 
3.75,  the  former  is  more  likely  to  be  correct 
ana  this  is  confirmed  by  the  calculations  for 
the  test  condition  with  values  of  K/O^/Ci  of 
0.80/54.5/1.39  (Fig  la).  Assumption  (a)  was 
therefore  adapted  for  all  subsequent  calcula¬ 
tions  for  this  cascade. 

Tne  choice  of  stream  tube  thickness  variation 
is  a  good  exanple  of  a  situation  where  there 
is  uncertainty  in  the  ope r ’.mental  data;  and 
hence  where  comparisons  against  other  nethods 
run  with  the  same  assumptions  would  be  parti¬ 
cularly  useful. 

An  89  x  16  calculating  grid  was  employed  for 
all  calculations  of  the  V2  cascade,  with  51  x 
16  points  equally  spaced  in  the  axial  and 
tangential  directions  within  the  blade  row. 
Outside  the  row,  the  axial  spacing  was 
increased  by  a  constant  factor  to  give  up¬ 
stream  and  downstream  boundaries  for  the 
calculation  domain  0.8  axial  chords  from  tire 
blade  row. 

The  twelve  VG18  test  conditions  can  be 
grouped  info  three  categories:  (a)  choked 
flow,  (b)  operation  near  optimum  incidence 
and  (c)  stalled  flow.  Tne  three  choked  flow 
conditions  are:-  0.85/49.5/1.20, 
0.80/47.5/1.17,  and  0.80/49.5/1.27.  They 
demonstrate  that,  relative  to  a  nominal 
unchoked  condition  of  0.80/49.5/1.20,  choked 
flow  can  be  produced  by  snail  changes  in  Mj, 
Oj,  or  n.  S1BYL2  successfully  predicts  this. 
For  example,  the  predicted  tech  number 
contours,  blade  surface  tech  nunber  distri¬ 
butions  and  development  of  the  boundary  layer 
displacement  thickness  and  shape  factor  at 
0.85/49.5/1.20  are  shown  in  Fig  lb,  together 
with  the  test  data.  There  is  good  agreement 
between  prediction  ard  test,  except  for  the 
levels  of  shape  factor  near  the  trailing  edge 
of  the  suction  surface. 

The  point  at  0.80/49. 571. 12  is  typical  of  the 
seven  operating  conditions  near  optimum 
incidence  with  measured  loss  coefficients 
below  5%.  The  predicted  tech  number 
contours,  blade  surface  tech  number  distri¬ 
butions,  and  boundary  layer  development  are 
shown  in  Fig  lc,  together  with  the  test  data. 
It  can  be  seen  that  the  overall  performance 
parameters  match  quite  well,  although  the 
discrepancy  of  1.3°  on  exit  flow  angle  is 
higher  than  usual.  There  is  also  close 
agreement  between  the  predicted  inviscid  tech 
numbers  on  the  blade  surfaces  and  the  values 
deduced  frcm  the  measured  surface  static 
pressures.  The  agreement  between  predi^ion 
arid  test  for  the  suction  surface  boundary 
layer  development  is  reasonable  qualitati¬ 
vely,  but  the  predicted  values  of  shape 
factor  are  higher  than  the  test  values  in  the 
separated  flow  region  near  the  trailing  edge, 
while  the  predicted  values  of  displacement 
thickness  near  mid  chord  are  lower  than  the 
test  values.  These  differences  could  be  due 
TM  P  1204 


to  many  reasons.  Seme  of  the  more  ob\iou3 
ones  are: 

i  errors  in  the  prediction  of  trans¬ 
ition  between  laminar  and  turbulent  flow 

ii  the  effect  of  liigh  free-stream 
turbulence  on  the  turbulent  boundary 
layer  -  this  is  not  included  in  the 
present  calculations 

iii  the  general  uncertainty  in  predic¬ 
ting  separated  turbulent  boundary  layers 

iv  errors  in  measuring  very  thin 
boundary  layers  -  the  measured  displace¬ 
ment  thickness  is  only  0.4  mm  at  42% 
chord 

v  errors  in  measuring  unsteady  flows 
with  a  pneumatic  probe 

vi  errors  in  measuring  reverse  flows 
with  a  forward  facing  probe 

The  level  of  agreement  achieved  in  this  case 
is  typical  of  that  for  the  test  conditions 
near  optimum  incidence. 

The  remaining  two  operating  conditions  - 
0.80/52.5/1.12  a»d  0.80/54.5/1.23  -  are 
stalled,  with  measured  loss  coefficients  of 
over  10% .  As  previously  noted,  such 
conditions  can  be  identified  by  the  fact  that 
maximum  blade  loading  occurs  at  the  leading 
edge.  The  user  can  then  specify  that  the 
suction  surface  boundary  layer  starts  as  a 
turbulent  layer,  with  a  momentum  thickness 
Reynolds  number  of  about  500.  A  solution 
with  the  boundary  layer  tripped  in  this  way 
for  the  condition  0.80/52.5/1.12  is  shown  in 
Fig  Id.  Tlie  results  for  surface  tech  number 
and  suction  surface  displacement  thickness 
agree  quite  well  but  the  predicted  shape 
factors  for  the  suction  surface  are  much 
higher  than  deduced  from  the  boundary  layer 
measurements . 

In  this  context  it  is  interesting  to  consider 
the  loss  produced  by  applying  a  compressible 
flow  mixing  calculation  to  the  measured 
suction  surface  boundary  layer  parameters  at 
98%  chord.  At  the  two  stalled  flow  condi¬ 
tions  this  loss  exceeds  the  total  loss 
deduced  from  the  wake  traverses  qy  about  25%. 
Assuming  that  the  suction  surface  loss  should 
be  about  the  same  proportion  of  the  total 
loss  as  at  tie  other  test  conditions  (ie  80%) 
implies  that  the  momentum  thickness  should  be 
about  half  of  the  stated  value.  This  would 
obviously  be  in  much  better  agreement  with 
the  predicted  values  of  momentum  thickness 
and  shape  factor. 

An  alternative  approach  to  tripping  the 
boundary  layer  would  be  to  employ  a  grid  with 
iruch  closer  axial  spacing  near  the  leading 
edge  and  to  assume  that  all  boundary  layers 
are  turbulent  This  would  give  solutions 
broadly  comparable  with  the  Navier-Stokes 
solutions  for  the  V2  cascade  shown  in  Ref  10. 
However,  a  grid  giving  better  modelling  of 
the  leading  edge  region  would  be  needed  to 
obtain  grid- independent  solutions.  Also  the 
process  of  laminar  separation,  transition  and 
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turbulent  reattachment  for  high  incidence 
flows  needs  to  be  understood  better  so  that 
more  realistic  models  of  the  viscous  flow 
(either  integral  boundary  layer  calculations 
or  turbulence  models)  can  be  developed. 

3.2  ARL  SL19  Transonic  Ccnpre33or  Cascade 

The  ARL  SL19  transonic  compressor  cascade  was 
derived  from  a  rotor  section  of  a  single 
stage  compressor  11  with  a  design  point 
pressure  ratio  of  1.9.  It  is  a  thin,  high 
stagger  (56.9°),  low  camber  (-1.9°)  section 
with  a  design  point  inlet  Mach  number  of 
1.616.  Cascades  with  similar  blades  have 
been  tested  in  facilities  at  Detroit  Diesel 
Allison  (DDA)  DFVLR  »  and  CNERA  >*,  and 
one  test  point  from  each  series  was  selected 
by  18 . 

A  guide  to  the  relevance  of  each  operating 
condition  can  be  obtained  by  computing  the 
pressure  ratio  and  efficiency  of  a  corres¬ 
ponding  rotor  section  from  the  measured 
overall  cascade  perfonran.ee  -  for  this 
purpose  the  rotor  section  was  assumed  to  have 
zero  inlet  absolute  flow  angle  and  constant 
stream  surface  radius.  The  results  for  all 
the  DDA  test  points  taken  at  the  design  inlet 


tech  number 

12  were  as  follows: 

:- 

Cascade 

Rotor 

Rotor 

static 

total 

adiabatic 

pressure 

pressure 

efficiency 

ratio 

ratio 

% 

1.220 

1.086 

67.8 

1.468 

1.180 

77.0 

1.672 

1.318 

81.8 

1.870 

1.489 

84.0 

2.036 

1.674 

84.8 

2.097 

1.727 

83.8 

2.220 

1.801 

81.1 

2.300 

1.868 

82.2 

This  clearly 

shows  that  the  points  with 

higher  static  pressure  ratios  are  of  greater 
practical  interest  to  the  compressor  engi¬ 
neer,  and  it  is  these  which  the  S1BYL2  method 
was  developed  to  predict.  Of  the  W318  test 
conditions  only  the  DFVIR  point  falls  into 
this  category. 

For  the  S1BYL2  calculation  of  the  DFVLR  test 
point  the  stream  tube  thickness  was  assured 
to  be  constant  upstream  of  the  cascade 
leading  edge  and  to  reduce  linearly  with 
axial  distance  to  the  measured  value  at  the 
traverse  plane  47%  of  gap  axially  downstream 
of  the  blade  trailing  edge.  A  119  x  21 
calculating  grid  was  employed  for  the  time 
marching  part  of  the  calculation,  with  81  x 
21  points  within  the  blade  row  equally  spaced 
in  both  axial  and  tangential  directions.  The 
test  value  of  inlet  tangential  velocity  was 
specified  and  the  operating  condition  for  the 
solution  was  chosen  such  that  the  main 
pressure  rise  on  the  suction  surface  occurred 
close  to  the  measured  position. 

The  predicted  tech  nuirber  contours,  blade 
surface  tech  number  distributions  and  bound¬ 
ary  layer  development  for  the  design  blade 
shape  are  shown  in  Fig  2,  together  with  the 
Mach  numbers  deduced  from  the  measured 


surface  static  pressures  and  the  predicted 
inviscid  total  pressures.  The  agreement 
between  prediction  and  test  is  generally 
good.  Note  that  the  calculation  indicates 
the  presence  of  strong  shock  wave/boundary 
layer  interactions  on  both  the  suction  and 
pressure  surfaces  of  the  blade.  It  is 
possible  that  the  pre-shock  boundary  layer  on 
the  pressure  surface  could  be  laminar, 
although  transition  has  been  forced  upstream 
of  the  shock  in  the  S1BYL2  solution,  and  this 
might  explain  the  discrepancies  between 
prediction  and  test  in  this  region.  The 
agreement  between  measured  and  predicted 
overall  performance  parameters  i3  reasonable 
-  the  predicted  unique  incidence  inlet  angle 
and  the  exit  flow  angle  are  both  1°  too  high, 
static  pressure  rise  agrees  to  about  5%,  and 
total  pressure  loss  i3  10%  compared  with  a 
measured  value  of  12%. 

The  predictions  are  quite  sensitive  to  the 
assumptions  made  about  stream  tube  thickness. 

For  example,  taking  a  constant  value  upstream 
of  the  initial  shock  wave  reduces  the  predic¬ 
ted  inlet  flow  angle  by  0.8°,  while  assuming 
a  linear  variation  between  the  plane  of  the 
upstream  static  pressure  tappings  and  the 
downstream  traverse  plane  increases  the  value 
by  0.8°.  There  i3  also  sene  uncertainty 
about  the  appropriate  blade  shape  to  use  for 
the  calculations,  since  the  two  instrumented 
blades  in  the  cascade  both  differed  slightly 
from  the  design  intent  and  seme  blade  deform¬ 
ation  due  to  the  aerodynamic  loads  wa3 
observed  during  testing  15 .  using  the  measu¬ 
red  shape  of  the  suction  surface  blade  given 
in  Ref  2  reduces  the  predicted  inlet  flow 
angle  by  0.7°,  whereas  the  shape  for  the 
pressure  surface  blade  is  nuch  closer  to  the 
design  intent. 

Cnee  again  calculations  with  a  number  of 
different  codes  would  be  instructive  and 
might  enable  a  consensus  view  to  be  reached 
of  the  appropriate  boundary  conditions  to  be 
applied. 

4,  HIGH  SPEED  OCMPREgSOR  TEST  CASES 

As  mentioned  above,  the  S1-S2  system 
incorporating  S1BYL2  has  been  applied  to  all 
three  high  speed  compressor  test  cases 
selected  by  WG18  viz: 

i  E/CO-2  NASA  Transonic  Rotor  67 

ii  E/QO-3  RR  HP9  Subsonic  Stage 

iii  E/CO-4  DFVIR  Transonic  Stage 

The  system  was  kept  the  same  throughout  and 
similar  calculating  grids  were  applied.  Nine 
SI  sections  ware  used  to  model  each  blade  row 
and  the  grids  used  for  the  time  marching  part 
of  the  S1BYL2  calculations  were  similar  to 
those  used  for  the  cascades  in  section  3 
above.  The  S2  solutions  typically  employed 
21  streamlines  and  there  were  4  interblade 
calculating  planes  between  the  leading  and 
trailing  edge  planes  of  each  blade  row. 

The  inlet  total  pressure  profiles  for  each 
compressor  were  based  closely  on  the  measured 
data,  except  that  no  attempt  was  made  to 
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achieve  zero  values  of  velocity  at  the  walls. 
Snail  amounts  of  annulus  wall  blockage,  about 
1%  for  the  two  transonic  fans  and  2 Vt%  for 
HP9,  were  included  to  allow  for  this.  The 
losses  used  to  nodal  end  effects  were  assumed 
to  affect  the  outer  regions  of  each  row:  the 
loss  coefficients  were  specified  to  drop  from 
15%  at  the  annulus  walls,  to  zero  at  a 
distance  of  1/3  rotor  chord  from  the  walls. 
This  assumption  is  in  line  with  that  made  for 
previous  S1-S2  calculations  of  the  rotor  row 
of  transonic  civil  fans6  “.  If  none  than 
one  stage  were  being  considered,  it  would 
obviously  oe  necessary  either  to  include  the 
effect  of  spanwise  mixing  within  the  S2  part 
of  the  calculation,  or  to  use  a  more  even 
distribution  of  extra  loss. 

4.1  NASA  Transonic  Rotor  67 

Rotor  67  is  the  fir3t  stage  rotor  of  a  two 
stage  fan.  It  was  tested  as  an  isolated  row 
to  provide  data  for  comparison  with  numerical 
predictions  which  were  free  frctn  circum¬ 
ferential  variations  due  to  stationary  blade 
rows.  The  rotor  has  an  aspect  ratio  of  1.6 
based  on  average  span/mid  height  chord,  and 
at  inlet  the  tip  diameter  is  514  mm  with  a 
hub/tip  ratio  of  0.375.  The  design  point 
conditions  are  a  tip  speed  of  429  m/s,  an 
inlet  tip  relative  Mach  number  of  1.38,  a 
rotor  ptessure  ratio  of  1.63  and  a  mass  flow 
rate  of  33.25  kg/s.  Laser  anemcrretry  was 
extensively  used  to  measure  details  of  the 
internal  flew  fields,  and  radial  traverses 
with  a  caibination  probe  were  carried  out  to 
determine  the  inlet  and  exit  variations  of 
total  and  static  pressure,  total  temperature 
and  flow  angle. 

Two  operating  conditions  were  selected  by 
W318  -  one  near  peak  efficiency  and  the  other 
near  stall.  The  flow  at  both  conditions  was 
predicted  using  the  S1-S2  system  and  same  of 
the  results  are  shown  in  Figs  3  and  4.  Con¬ 
sidering  firstly  the  overall  performance,  the 
mass  flow  rate  at  the  pressure  ratio  for  the 
peak  efficiency  condition  was  predicted  to  be 
34.0  kg/s,  about  l/:%  below  the  test  value. 
Predicted  adiabatic  efficiency  was  92.3%, 
compared  with  a  test  value  of  93%.  At  the 
mass  flow  for  the  near  stall  point,  the 
predicted  pressure  ratio  and  efficiency  were 

1.727  and  91.6%,  compared  with  test  values  of 

1.728  and  90.1%.  Tnere  is  also  good  agree¬ 
ment  between  the  predicted  and  measured 
radial  profiles  of  total  pressure  and 
temperature  at  rotor  exit,  as  shown  in  Fig  3. 
(The  values  plotted  are  relative  to  the  inlet 
conditions  at  mid  span.) 

The  predicted  Mach  number  contours  at  the  two 
conditions  for  the  blade-to-blads  sections  at 
30%,  70%  and  90%  span  from  the  hub  are  shown 
in  Fig  4,  together  with  the  predicted  blade 
surface  Mach  numbers  and  boundary  layer 
development.  The  contours  may  be  compared 
directly  with  the  cot responding  plots  deduced 
frcrn  the  test  data  given  in  Fig  17  of  Pef  2. 
The  agreement  is  quite  good  for  the  outer 
sections,  with  the  S1-S2  system  correctly 
predicting  the  dual  shock  system  noted  in  Ref 
l7  at  the  peak  efficiency  condition,  and  the 
single  normal  shock  wave  just  detached  frem 
the  blading  at  the  near  stall  condition. 


However,  the  calculation  predicts  a  more 
distinct  shock  in  the  uncovered  passage  at 
30%  blade  height  for  both  conditions  than  was 
measured  with  the  laser  anemometer. 

The  predicted  boundary  layer  development  for 
the  3uction  surface  indicates  a  trailing  edge 
separation  for  the  sections  over  the  inner 
third  of  the  3pan,  with  the  shock  wave/ 
boundary  layer  interaction  near  mid  chord 
becoming  significant  and  causing  sane  separ¬ 
ation  over  the  outer  half  of  the  blade.  The 
interaction  is  sufficiently  strong  at  the 
peak  efficiency  condition  to  cause  carplete 
separation  of  the  boundary  layer  downstream 
of  the  shock  for  the  tip  sections. 

4.2  RR  HP9  Subsonic  Compressor  Stac.e 

The  Rolls-Royce  .'P9  rig  is  a  single  stage 
model  of  one  of  t  e  later  stages  in  a  civil 
HP  compressor,  it  therefore  has  a  much 
higher  inlet  hub/tip  ratio  (0.84)  and  lower 
rotor  aspect  ratio  (0.9)  than  the  other  two 
W318  compressor  cases.  Inlet  spoiler  rings 
are  employed  to  generate  a  suitable  inlet 
velocity  profile.  The  type  "RD"  case 
selected  by  W318  has  conventional  blade 
profiles,  a  design  stage  pressure  ratio  of 
1.24,  a  mass  flow  rate  of  9.1  kg/s  and  a  tip 
speed  of  251.3  m/s  at  ISA  inlet  conditions. 
Tip  diameter  is  518.2  ran. 

All  three  operating  conditions  -  maximum 
flow,  mid  chic  and  near  surge  -  were  studied 
using  the  S1-S2  system.  The  predictions  at 
maximum  flow  and  mid  chic  were  straight¬ 
forward,  but  the  initial  attempts  for  the 
near  3urge  point  failed.  The  high  incidences 
for  the  SI  sections  nearest  the  hub  caused 
very  high  levels  of  profile  loss  to  be 
predicted  for  this  region:  these  resulted  in 
even  more  arduous  flow  conditions  being 
calculated  for  the  next  S1-S2  loop  by  the  S2 
part  of  the  system  and  eventually  the  S1BYL2 
calculations  failed  to  converge.  The  present 
S1-S2  technique  is  probably  over-3ensitive  to 
local  flow  problems  like  this,  since  there  is 
no  al'^w  nee  for  spanwise  movement  of  low 
energy  ooundary  layer  flows  -  inclusion  of 
spanwise  mixing  in  the  S2  calculation  should 
help  in  this  respect.  Also,  the  assumption 
that  stream  surfaces  remain  axisynuetric 
means  that  the  relief  of  tire  endwall  flow 
conditions  due  to  three-dimensional  flows  “ 
is  neglected,  a  fundamental  limitation  of  the 
axisyrmetric  S1-S2  approach.  To  enable  a 
solution  to  be  obtained  for  HP9  at  the  near 
surge  mass  flow,  SI  calculations  were  carried 
out  only  between  10%  and  90%  blade  heights 
(the  other  S1-S2  solutions  included 
calculations  for  sections  within  5%  height  of 
the  walls)  with  extrapolation  used  to  obtain 
the  conditions  for  stream  surfaces  closer  to 
the  walls. 

Predicted  Mach  number  contours  and  distribu¬ 
tions  and  boundary  layer  development  for  the 
rotor  and  3tator  mud  height  sections  at  the 
three  operating  conditions  are  shown  in  Fig 
5.  These  demonstrate  the  range  of  operating 
incidences  which  are  to  be  expected  at  diffe¬ 
rent  points  along  a  constant  speed  character¬ 
istic  for  a  single  stage  compressor  -  high 
positive  blade  loadings  at  the  leading  edge 
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near  surge,  roughly  zero  loa  'lug  at  mid  chic 
and  high  negative  loadings  at  i.vaxinum  flow. 
Th'-!  predicted  boundary  layer  behaviour  is 
typical  of  blade  sections  with  circular  arc 
camber  lines.  On  the  suction  surface  the 
shape  factors  are  relatively  low  over  most  of 
the  chord,  indicating  attached  boundary 
layers.  However,  there  is  rapid  growth  near 
the  trailing  edge  where  the  boundary  layer 
has  become  too  thick  to  withstand  t)>e 
constant  diffusion  associated  with  the 
circular  arc  canter  line  and  separates .  On 
the  pressure  surface  the  boundary  layer 
remains  attached  throughout,  apart  from  a 
possible  laminar  bubble  near  the  leading  edge 
at  the  maximum  flow  point.  No  interblade 
measurements  are  available  for  comparison 
with  these  data. 

There  is  3cme  uncertainty  about  the  rotor 
exit  conditions  for  this  case  !,  and  so  the 
comparisions  shown  between  prediction  and 
te3t  (Fig  6)  are  for  the  circumferentially 
averaged  conditions  downstream  of  the  stator. 
The  main  discrepancies  are  that:- 

i  the  temperature  rises  near  the 
outer  casing  are  overestimated  -  this 
could  be  due  to  the  relatively  crude 
allowance  for  end  wall  effects,  or  to 
neglecting  spanwise  mixing  effects  in 
the  S2  throughflow  calculation. 

ii  the  predicted  efficiencies  at  the 
near  surge  condition  are  much  higher 
than  measured  -  this  is  probably 
associated  with  the  problems  in 
predicting  high  incidence  flows  noted  in 
section  3.1  above.  Improved  agreement 
could  be  obtained  by  "tripping"  the 
suction  surface  boundary  layers  at  the 
blade  leading  edge. 

Apart  frcm  this  the  agreement  obtained  is 
reasonable,  particularly  when  the  uncertainty 
due  to  the  manufacturing  tolerance  of  0.75° 
on  blade  stagger  angle  is  considered.  If  the 
prediction  at  the  maximum  flow  point  is 
repeated  with  the  rotor  and  stator  rows 
restaggered  open  by  this  amount,  then  at  the 
same  exit  flow  function  there  are  increases 
of  0.9%  on  mass  flow,  0.012  on  pressure  ratio 
ana  1.2%  on  efficiency. 

4.3  DF7IK  Transonic  Compressor  Stage 

The  DFVLR  transonic  fan  has  an  inlet  tip 
diameter  of  398  mm,  an  inlet  hub/tip  ratio  of 
0.5  and  a  rotor  aspect  ratio  of  1.7.  The 
design  point  pressure  ratio  is  1.51  at  a  mass 
flow  rate  of  17.3  kg/s  and  a  blade  tip  speed 
of  421  m/s.  The  rotor  has  been  tested  with 
two  different  stator  rows:  the  laser  aremo- 
metry  measurements  of  the  rotor  flow  field 
were  taken  with  the  NACA  65/60  stator  row  and 
this  is  the  configuration  which  has  been 
studied  using  the  S1-S2  system. 

Predictions  were  carried  out  at  the  same  exit 
mass  flow  functions  as  the  two  test  condi¬ 
tions  selected  by  W318  and  the  predicted  and 
measured  radial  distributions  of  performance 
are  compared  in  Fig  7.  The  agreement  is  much 
poorer  than  for  the  other  two  compressors 
The  main  problem  seems  to  be  that  the  S1-S2 


system  over-estimates  the  rotor  efficiency  by 
about  4%  at  both  conditions,  and  this  results 
in  significantly  higher  pressure  ratios  being 
oredicted  at  a  given  mass  flow.  The  magni¬ 
tude  of  this  discrepancy  is  disappointing, 
gi\  in  the  level  of  agreement  on  overall 
performance  for  the  other  two  cases  and  for 
otlw-r  transonic  fans.  In  particular  the 
S1-S2  system  has  not  identified  any  reason 
tne  peak  efficiency  of  the  DFVIR  rotor 
bei’ig  over  5%  lower  than  for  NASA  rotor  67, 
which  has  a  similar  pressure  ratio,  tip  speed 
aspect  ratio.  It  will  be  interesting  to 
see  whether  other  CFD  methods  can  explain 
chis  difference. 

Despite  the  poor  agreement  on  overall  per¬ 
formance,  the  predicted  Mach  number  contours 
for  the  rotor  at  18%,  45%,  68%  and  89%  blade 
heights  are  in  close  agreement  with  those 
given  in  Figs  10  and  li  of  Ref  2  (See  Fig  8) . 

At  peak  efficiency  the  solution  correctly 
predicts  an  oblique  leading  edge  shock  near 
the  rotor  tip  with  a  maximum  suction  surface 
Mach  number  of  about  1.4.  At  68%  height  the 
shock  is  becoming  more  normal,  and  by  45% 
height  it  is  detached  frcm  the  blacte  leading 
edge.  The  predicted  diffusion  of  the  flow 
within  the  blade  passage  also  matches  the 
test  data  well.  The  predicted  suction 
surface  boundary  layer  development  shows  a 
gradual  progression  with  blade  height, 
roughly  similar  to  that  for  NASA  rotor  67. 

Near  the  hub,  separation  occurs  only  near  the 
trailing  edge;  by  about  45%  height  the  shock 
wave/ooundary  layer  interaction  has  become 
strong  enough  to  cause  a  small  separation 
bubble  at  about  mid  chord,  and  the  extent  of 
the  separation  at  the  trailing  edge  has 
reduced;  at  68%  height  the  mid  chord  separa¬ 
tion  has  increased  and  the  trailing  edge 
separation  has  disappeared;  and  at  89%  the 
shock  wave/boundary  layer  interaction  is 
strong  enough  to  cause  complete  separation 
for  the  rear  part  of  the  blade. 

At  the  point  nearer  surge,  the  shock  wave  has 
moved  forward  and  it  is  just  detached  frcm 
the  leading  edge  at  89%  height  in  both 
predicted  arid  measured  results.  The  shock 
becomes  further  detacher  at  68%  and  45% 
heights,  anu  at  18%  there  is  only  a  small 
region  near  the  leading  edge  where  the  flow 
is  supersonic.  The  predicted  boundary  layer 
behaviour  is  similar  to  that  at  peak  effici¬ 
ency  except  that  Ute  shock/boundary  layer 
interaction  near  the  tip  is  predicted  to  be 
slightly  weaker  and  the  boundary  layer 
reattaches  before  the  trailing  edge.  How¬ 
ever,  given  the  high  effective  incidences 
indicated  over  the  inner  half  of  the  blade, 
it  is  likely  that  the  actual  suction  surface 
boundary  layers  are  rather  thicker. 

5.  CCTCTUDIN3  REMARKS 

The  S1BYL2  inviscid-vincous  interaction 
method  has  now  been  applied  to  all  the  high 
speed  compressor  cascade  test  cases  selected 
by  AGARD  W318,  and,  linked  with  a  streamline 
curvature  throughflow  method  in  an  S1-S2 
system,  to  all  the  high  speed  compressor 
cases.  The  predictions  detained  are 
generally  in  encouraging  agreement  with  the 
experimental  data,  both  in  terms  of  the 
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features  of  the  internal  flows  and  the 
overall  performance  parameters . 

For  the  ccmpresscr  cascades  the  blade 
choking/unique  incidence  angles  and  deviation 
angles  are  predicted  to  about  ±  1°,  the 
static  pressure  rise  to  ±  4%  and  (less 
satisfactorily)  the  total  loss  coefficient  to 
±  2%  of  inlet  dynamic  lead,  for  operation 
below  stalling  incidence.  The  increased 
losses  wiiich  occur  at  high  positive  incidence 
are  not  fully  predicted  automatically,  but 
such  situations  can  be  detected  by  examining 
the  effective  incidence  predicted  for  the 
blade.  It  is  believed  that  improvements  are 
needed  both  in  the  viscous  modelling  and  in 
the  numerical  resolution  of  the  leading  edge 
region  before  this  type  of  flow  can  be 
predicted  correctly. 

For  the  compressor  test  cases,  the  predicted 
results  for  the  NASA  rotor  67  and  for  the  KP9 
stage  agree  well  with  the  measured  results. 

In  particular,  peak  efficiency  levels  are 
estimated  to  within  about  1%.  The  results 
for  the  DFVUR  stage  are  less  satisfactory: 
the  rotor  efficiency  is  over-estimated  by 
about  4%,  despite  close  agreement  between  the 
measured  and  predicted  blade-to-blade  Mach 
numbers.  Thus  the  present  method  fails  to 
explain  the  difference  of  over  5%  in  measured 
rotor  efficiency  betwsen  the  NASA  67  and 
DFVLR  rotors.  It  will  be  interesting  to  see 
whether  other  methods  are  more  successful. 

The  present  applications  of  the  RAE  S1BYL2 
code  to  the  W318  compressor  test  cases  have 
been  valuable,  both  in  evaluating  the  code 
and  in  gaining  further  understanding  of  the 
test  cases.  The  author  would  like  to 
congratulate  the  Working  Group  for  taking  the 
lead  in  identifying  this  set  of  cases,  and 
he  would  encourage  others  to  carry  out  and 
publish  corresponding  predictions. 
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Fig  2  S1BYL2  prediction  For  ARL  Sl’9  cascade  (symbols  indicate  test  data) 


a.  Padial  profiles  near  peak  efficiency  b.  Radial  profiles  near  stall 


Fig  3  S1-S2  predictions  For  NASA  transonic  rotor  67 
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Fig  G  S1-S2  predictions  For  Rolls-Royce  HP9  stage 


a.  Radial  profiles  near  peak  efficiency 


Fig  7  S1-S2  predictions  For  DFVLR  transonic  stage 
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